Chiral rf-wave superconductor with nonzero center-of-mass pair momentum 
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Hypothetical topologically nontrivial superconducting state of two-dimensional electron system 
is discussed in connection with the problem of pairing with large center-of-mass pair momentum 
under predominant repulsive screened Coulomb interaction. Direct numerical solution of the self- 
consistency equation exhibits two nearly degenerate order parameters which can be formally referred 
to d x s 2 and d xy orbital symmetry. Spontaneous breaking of the time-reversal symmetry can mix 
these states and form fully gapped chiral d + id superconducting state. 
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In recent years, condensed matter physics has 
significantly focused on studies of peculiar states of 
matter, such as two-dimensional (2D) topological 
insulators and superconductors i Both time-reversal 
(TR) invariant and TR breaking topological 
superconductors have attracted a lot of interest, in 
particular, because of their potential applications. 
TR breaking superconductors are classified by an 
integer topological invariants 2 similar to those used 
for classification of quantum Hall states^ The simplest 
chiral topological triplet superconductor with p x + ip y 
orbital symmetry was considered by Read and Green 4 
and predicted to exist in S^RuCU by Mackenzie and 
Maeno£ 

Degeneration of the d , 2 and eL„, ordered states 
inherent in doped graphene monolayer has recently 
considered by Nandkishore et al£ as possible origin of 
a rise of a singlet chiral superconducting (SC) state with 
d + id orbital symmetry. Such complex order parameter 
was suggested by Laughlin 7 to connect the TR broken 
symmetry and the low-temperature phase transition 
observed in Bi2Sr2CaCu20s in external magnetic field. 8 
In such a case, the d xy component of the SC order 
parameter turns out to be field-induced. Similar phase 
transition was observed in Ni-doped Bi2Sr2CaCu20*8 in 
zero external magnetic fields Balatsky 10 pointed out 
that, in the presence of magnetic impurities, the d x2 _ y2 
superconductor can exhibit a transition exactly into 
the d + id state due to a coupling between impurity 
magnetization and the d xy component of the order 
parameter. 

In this Communication, we report a possibility of a rise 
of d + id chiral state in high-temperature SC cuprates. 

Recently introduced concept of SC pairing with large 
pair momentum under screened Coulomb repulsion 11 
offers an explanation of principal features of cupratc 
superconductors: 1) a checkerboard real-space ordering 
observable in the SC stated 2 - can be directly related to 
the pair momentum, comparable with reciprocal lattice 
spacing; 13 2) pseudogap (PC) stated with a broad region 
of SC fluctuations above the transition temperature T c 
can be explained by a rise of quasi-stationary states of 
pairs with large momentum due to real-space oscillations 
of the screened Coulomb potential^ 3) high-energy 



effects observable in optical experiments^ can be related 
to the electron-hole asymmetry that becomes apparent 
in the SC state of the cuprates. 

The order parameter originating from SC pairing with 
large momentum is nonvanishing in the interior of a part 
of the Brillouin zone (domain of kinematic constraint) 
due to the fact that, at T = 0, the momenta of both 
particles composing a pair should be either inside or 
outside the Fermi contour (FC). This order parameter 
turns out to be appreciably nonzero inside vicinities of 
nested segments of the FCj^ 5 . In the case of SC cuprates, 
such segments correspond to antinodal region of the 
Brillouin zone. 12 

High values of T c and specific isotope effect manifested 
in the cuprates^ 7 - show that, together with the 
repulsive Coulomb pairing interaction, one should take 
into account the attractive contribution owing to 
electron-phonon interaction (EPI) , including the forward 
scattering effect ^ 

In the antinodal region, phonon assisted Coulomb 
pairing with large momentum can predominate 
over conventional phonon-induced pairing with zero 
momentum that prevails only in the nodal region. Then, 
superconductivity at low temperatures should exist as 
a biordered state formed by the condensates of pairs 
with large and zero momenta in the antinodal and 
nodal regions of the momentum space, respectively. 15 
On the contrary, the SC order just below T c should 
be determined by the pairing with large momentum. 
Therefore, such an order should arise only in the 
antinodal region. 

Considerable enhancement of T c , observable in the 
cuprates, can be qualitatively related to specific phonon 
induced "symmetrization" of real-space oscillations of the 
screened Coulomb potential* 1 ^ The SC order parameter 
A(fe; K), where k and K are relative motion and center- 
of-mass momenta of SC pair, respectively, should be 
obtained as a self-consistent solution to the mean-field 
gap equation with the momentum representation of such 
a symmetrized potential U{k,k';K). The gap equation 
at T — can be written as 
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Here, 2e + (k; K) = e(k + ) + e(k_) is kinetic energy of 
a singlet pair composed of the particles with momenta 
k ± = K/2±k, the summation is taken over momenta k' 
belonging to the domain of kinematic constraint relevant 
to given pair momentum K. 

In the case of pairing with K ^ 0, in common 
with the well-known Fulde-Ferrel-Larkin-Ovchinnikov 
(FFLO) problem ; 20 ! 21 the order parameter can be 
represented in the form of either a running wave,— A ~ 
expiKR, or a standing wave,— A ~ cos (K R), that is 
as a symmetric superposition of running waves with pair 
momenta ±_KT. Here, R is center-of-mass radius- vector of 
the pair. It should be noted that, unlike the FFLO state, 
the SC state with nonzero center-of-mass momentum 
considered here arises without external magnetic field 
and, generally speaking, preserves TR symmetry. 

One can expect that, along with the symmetric 
superposition, antisymmetric superposition of the same 
waves, A ~ sin (KR), could be a solution to the gap 
equation as well. Both symmetric and antisymmetric 
solutions are defined in common domain of kinematic 
constraint which has to be constructed as the union of 
the domains for running waves with ±K. 

Because of the crystal symmetry of the system, FFLO 
order parameter can be defined as a more complicated 
linear combination of running waves with equivalent 
momenta^ In a similar way, one can define zeroth-order 
approximation of the order parameter arising as a result 
of SC pairing with large pair momentum. 

In the case of the cuprates, tetragonal symmetry 
of Cu02 plane results in four crystal equivalent pair 
momenta: ±K and ±K' where K' is perpendicular 
to ±K . It is convenient to form standing waves as 
symmetric and antisymmetric (with respect to in-plane 
reflection from a line perpendicular to pair momentum) 
superpositions for each of two running waves with 
momenta, ±K and ±K', respectively. Then, the order 
parameter in the whole of the Brillouin zone can be 
written as a linear combination of these standing waves. 
Coefficients in such linear combinations specifies the 
orbital symmetry of the order parameters- 
Coefficients of like signs correspond to extended s- 
wave symmetry (the order parameter is invariant with 
respect to rotation by ir/2 about C 4 axis). In the case of 
coefficients of unlike signs, the order parameter reverses 
sign under rotation by it/ 2 and therefore can be referred 
to d-wave orbital symmetry. 

One can represent the order parameter by any of four 
linear combinations (two s-wave and two <i-wave) directly 
following from the gap equation. The SC ground state of 
the system should be expressed by the linear combination 
which has a lower free energy. 

For solving the gap equation, we present interaction 
energy U(k,k';K) as a sum of screened Coulomb 
repulsion, U s (k,k'), defined in the whole of the domain 
of kinematic constraint and EPI induced effective 
attraction which is assumed nonzero inside a narrow 
region enveloping the FC^ Width of this region in the 



momentum space is of the order of 2lo d /v f where uj d 
and v F are characteristic Debye frequency and Fermi 
velocity normal to the FC, respectively. We assume 
that attractive contribution into interaction energy is 
nonzero if and only if momenta of particles before and 
after scattering (k and k', respectively) both belong to 
this region. Also, we assume that this contribution is 
independent of momenta inside the region. 

Thus, U{k,k';K) = U s (k,k') - V, if both momenta 
belong to the region, U(k, k': K) = U s (k, k') when even 
if one of the momenta k and k' belonging to the domain 
of kinematic constraint does not belong to the region. 
One can assume that U(k, k'; K) = when even if one of 
the momenta does not belong to the domain of kinematic 
constraint. 

To solve the gap equation numerically, we approximate 
the interaction energy by a simple step function of n — 
k — fc'. 13 Step length in the direction of K is limited 
by the length of the nested segment of the FC. In the 
direction perpendicular to K , left step at k 1 corresponds 
to phonon-mediated decrease in energy taking into 
account the forward scattering effect 18 whereas right step 
at K r reflects the fact that values of the order parameter 
turn out to be very small at all points distant from 
the FC. Therefore, the order parameter turns out to be 
weakly sensitive to 

We use the electron dispersion that conforms to the 
FC observable in hole doped cuprates, 

s(k x , k y ) = 2t — 2t(cosk x + cosk y ) — At' cosk x cosk y — 
- 2t"{cos2k x + cos2k y ), (2) 

where t = 1 eV, t = 0.5 eV, t'/t = -0.3, t"/t = 0.14, and 
dimensionless components of momentum k i (i = x, y) 
vary within — tt < k i < n. 

Recently obtained numerical solution to the gap 
equation 13 corresponding to a symmetric superposition 
of two running waves is presented in Fig. 1 (only left half 
of the Brillouin zone is shown) in which the momenta 
of the running waves ±K are chosen as directed along 
k y axis. The order parameter A s (k; K) is characterized 
by intrinsic system of zero lines intersecting the FC. 
One can see that A s (fc; K) possesses distinct values only 
in a vicinity of nested segments of the FC. Similarly, 
one can obtain the symmetric superposition A s (k;K') 
corresponding to the running waves with ±K' . 

The symmetric superposition A s (k;K) can describe 
a stripe structure of the SC stated in particular, an 
emergence of the SC order parameter appearing in the 
PG state of the cuprates^ Recently, Berg et ali2i 
discussed similar striped SC state as a unidirectional 
pair-density wave (PDW) phase with periodic real-space 
dependence of the order parameter on the center-of-mass 
position. The coupling between such a PDW and other 
ordered states was also considered in the framework of 
Ginzburg-Landau theory^ 

To obtain the order parameter in the whole of 
the Brillouin zone, one should compose either s-wave 
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Phc. 1. Symmetric superposition A s (k; K) of running waves 
with wave vectors ±K directed along k y axis. Left panel: 
dependence of A s (k; K) on k y at k x — —0.7837T that 
corresponds to the position of the left vertical nested segment 
of the FC. Right panel: topology of A s (fc; K); only left half of 
the Brillouin zone is shown. White lines present the intrinsic 
zero lines on which A s (fe; K) = 0. White background between 
the Fermi contours shifted by ±K /2 with respect to initial 
position of the FC reflects the kinematic constraint. 
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Phc. 2. Antisymmetric superposition A a (k; K) of running 
waves with wave vectors ±_K" directed along k y axis. Left 
panel: dependence of A a (k; K) on k y at k x — — 0.7837T that 
corresponds to the position of the left vertical nested segment 
of the FC. Right panel: topology of A a (fc; K); only left half of 
the Brillouin zone is shown. White lines present the intrinsic 
zero lines on which A a (fc; K) = 0. White background between 
the Fermi contours shifted by ±K/2 with respect to initial 
position of the FC reflects the kinematic constraint. 



or d-wave linear combination of obtained symmetric 
superpositions,— 

A^(k)~A s (k;K)±A s (k;K'). (3) 

The first of them (A^" 1 ), corresponding to extended s- 
wave orbital symmetry, displays the intrinsic zero lines 
only, whereas the second one (Ai - -*), corresponding to 
d x 2_ y2 orbital symmetry, besides the intrinsic zero lines, 
displays four straight zero lines (nodal lines) along the 
diagonals of the Brillouin zone. 

In the cuprates, it seems that such orbital d-wave nodal 
lines are consistent with available experimental facts 
including angle-resolved photoemission spectroscopy 
(ARPES) data*2£. In particular, it is very likely that SC 
gap near the diagonals takes small values and can even 
vanish. On the contrary, intrinsic zero lines, situated close 
to the FC, hardly ever can be detected directly from 
ARPES measurements but they undoubtedly should 
become apparent in thermodynamical properties. 

In Fig. 2, we present new numerical solution to the 
gap equation, namely, the order parameter in the form of 
antisymmetric superposition of two running waves with 
opposite momenta, A a (k:K). 

This solution, just as symmetric superposition 
A s (k;K), has its own system of intrinsic zero lines. As 
one can see from Fig. 2, one of such lines turns out to be 
parallel to one of the coordinate axes. 

Extreme values of both A s (k;K) and A a (k; K) 
are concentrated in the antinodal regions near nested 



segments of the FC. It should be noted especially 
that, as results from direct numerical solution to 
the gap equation, both symmetric and antisymmetric 
superpositions exhibit comparable extreme values. 

The order parameter in the whole of the Brillouin zone 
should be presented as either s-wave or d-wave linear 
combination of antisymmetric superpositions, 

A^(k)~A a (k;K)±A a (k;K'). (4) 

The s-wave combination corresponding to plus sign in 
Eq. (HJ), besides some closed intrinsic zero lines, displays 
eight straight zero lines parallel both sides and diagonals 
of the Brillouin zone. Formally, such an order parameter 
can be referred to the so-called g-w&ve orbital symmetry 
discussed by Zhao»2& 

The d-wave combination corresponding to minus sign 
in Eq. (JJ), besides some closed intrinsic zero lines, 
displays four intrinsic straight zero lines directed along 
the sides of the Brillouin zone. These lines can be formally 
considered as the nodes of the order parameter with d xy 
orbital symmetry. 

It is the pairing interaction that makes the choice in 
favor of a certain set of coefficients in Exs. @ and ((3]). 
For example, near a spin-density-wave instability typical 
of the cuprate superconductors, paramagnon exchange 
induced component of the pairing interaction, sensitive 
to band structure and band filling, gives rise to singlet d- 
wave pairing with zero center-of-mass momentum. Such 
an interaction favors (suppresses) d x2 _ 2 (d xy ) channel. 27 
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As follows from Figures 1 and 2, the amplitudes of two 
superpositions, A s (fe; K) and A a (fc; K), are comparable. 
Therefore, under certain conditions, the order parameter 
with spontaneously broken TR symmetry, 

A( c )(fc)~A(-)(fc)+zA(-)(fc), (5) 

can be expected as a chiral fully gapped ground state of 
the system. The components of A^-'(fe), phased by n/2 
with respect to each other, should be formally referred to 
d-wave orbital symmetry, d x2 _ y2 and d xy , respectively. 

It should be noted that, in the problem of 
triplet pairing with large center-of-mass momentum, 
superpositions A a (k;K') and A a (k; K) themselves, 
can be considered as the p-wave order parameters 
corresponding to p x and p y symmetry, respectively. These 
superpositions can be used to form chiral triplet SC state 
with p x + ip y orbital symmetry. 

Values of the d xy component, Aq ~\k), should be 
nonzero for momenta belonging to the nodal region 
because of the proximity effect in the momentum space 
and a contribution of SC pairing with zero momentum 
into biordered SC stated (both ignored in Figures 1 and 
2). As a result, a finite SC gap should appear at the points 
that correspond to the nodes of the pure d^„ 2 state. 

x y 

The d xy order parameter is often accepted in a 
simple form, sinfc^sinfc^, that reveals maxima just on 
the diagonals. On the contrary, d + id order §5§ turns 



out be concentrated in antinodal vicinities of the nested 
segments of the FC (see Fig. 2). Therefore, the smallness 
of the d xy component of Ex. ([5]) in the nodal region 
can make it difficult to detect such a gap (for example, 
using ARPES technique). It should be noted, however, 
that recent ARPES data^ can be considered as an 
unambiguous evidence in favor to d 2 . 2 _ l ,2-wave-hke SC 
gap in optimally doped YBa2Cu307_a that exhibits a 
nonzero minimum of about 12meV along the nodal 
direction. 

Chiral SC ground state with the order parameter 
([5]) turns out to be topologically nontrivial. Indeed, by 
continuous deformation of the parameters of the mean- 
field Hamiltonian without opening a gap, this state 
can be transformed into that considered by Volovik: 29 
both states are topologically equivalent and should be 
characterized by the topological invariant J\f = ±2. 
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